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ABSTW^CT 

It has been suggested that surface waves nay be able to heat the solar 
corona. These waves can propagate into the corona and supply the tequired 
energies, and because they are linearly cooqpressive they can be dissipated 
by ion viscosity and electron heat conduct! m. In this paper we evaluate 
the damping of surface waves by viscosity and heat conduction. For the 
solar corona, it is found that surface waves dissipate efficiently only if 
their periods are shorter than a few tens of seconds and only if the back- 
ground magnetic field is less than about 10 Gauss. Heating of quiet coronal 
regions is possible if the coronal, waves have short periods, but they cannot 
heat regions of strong magnetic field, such as coronal active region loops. 


I. Introduction 


The mechanisms which convert the kinetic energy of the solar photosphere 
and convection zone into the thermal energy of the corona have been under 
investigation for some four decades. In spite of intensive and imaginative 
effort, no successful theoretical solution to the coronal heating problem 
has been offered (see recent reviews by Chiuderi , 1979; Hollweg , 1981a; 

Kuperus, lonson and Spicer , 1981; Wentzel, 1981; Withbr je, 1981). 

Coronal heating by MHD waves has been widely investigated. The MHD 
slow mode can probably be Ignored in this context, since it propagates too 
slowly to carry the required energy flux into the corona, subject to the 
constraints imposed by the observed amplitudes of non-thermal motions in 
the corona and underlying chromsphere (e.g. Athay and White , 1979a, b; Bruner , 
1978; Bruner and Poletto , 1981). 

The MHD Alfven mode has been investigated in this context for many 
years (e.g. Alfven , 1947; Osterbrock , 1961; Piddlngton , 1956). Even though 
\t is strongly reflected by the chromosphere corona transition region, the 
Alfven mode does seem capable of carrying the required energy flux into the 
corona from below, at least if the wave frequency is not too low (e.g. 

Hollweg , 1981b; Hollweg, Jackson, and Galloway , 1982; Leer, Holzer and Fla , 
1982; Zugzda and Locans , 1982). The constraints are more severe in active 
regions, but Hollweg (1981) and Zugzda and Locans (1982) have shown that 
energy can enter coronal active region loops if appropriate resonances in 
the coronal active region loops are excited. The effect is similar to anti- 
reflectance coatings on camera lenses; lonson ( ) has discussed it in terms 


of resonant LRC circuits. The difficulty with the Alfven mode is the 
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dissipation mechanism. Alfven waves are noncompressive (to lowest order), 
and they therefore do not couple to the radiation field, nor are they 
damped to any significant extent by viscosity or heat conduction. Alfven 
wpves are expected to be only weakly nonlinear in the corona, and nonlinear 
mechanisms are probably inadequate (e.g. Chin and Wentzel , 1972; Uchida and 
Kaburaki , 1974). A possible exception to this statement has been pointed 
out by Hollweg, Jackson and Galloway (1982), who noted that Alfven waves 
with periods shorter than a few minutes can steepen into shocks in the 
chromosphere, which subsequently enter the corona from below. Hollweg (1982a) 
has shown that the shocks can in principle dissipate rapidly enough to heat 
coronal holes or quiet coronal regions, but the mechanism fails on coronal 

m 

active region loops. 

The MHD fast mode has a number of appealing features. Its group velocity 
is large and it is therefore in principle capable of carrying a large energy 
flux in the corona. It is intrinsically compressive, and therefore subject 
to dissipation by viscosity, heat conduction and radiation, or by Landau and 
transit-time damping in the high-frequency limit where Coulomb collisions 
are ineffective. Habbal, Leer, and Holzer (1979) have pointed out that re- 
fraction can focus the wave energy into selected sites in the corona, thus 
offering a natural explanation of the observed "structuring'* of the corona. 
However, Habbal, Leer and Holzer simply postulate the existence of an appro- 
priate flux of fast waves at the coronal base. This approach may not be 
justified. Ostei brock (1961) has pointed out that fast waves should undergo 


severe refraction in the chromosphere, where the Alfven speed increases 


5 


rapidly with halght; the refraction turns the energy flux away from the 
corona. In a similar vein, Hollweg (1978) has pointed out that the 
known solar motions should lead primarily to fast waves which are 
evanescent In the corona. These conclusions have recently been put on 
firmer ground by Leroy and Schwarts (1982) and Schwarts and Leroy (1982) , 
who have Investigated the behavior of fast waves In a model solar atmos- 
phere which Includes the chromosphere, transition region, and corona; 
they conclude that fast waves can not supply the required energy flux Into 
the corona. 

A remaining possibility takes advantage of observations that the 
corona Is highly structured. It Is possible that the structuring occurs 

m 

In the form of very thin tangential discontinuities. Such structures would 
have thicknesses of only a few proton gyro-radll, which would be totally 
unresolved In any optical observations of the corona. However, It Is 
worth noting that such tangential discontinuities are very abundant In the 
solar wind (e.g. Burlaga , 1971; Neubauer and Barnstorf , 1981), occurring 
as frequently as 1 hour~^; If the solar wind can be used as a guiue, the 
corona may contain tangential discontinuities (TD's). It has been pointed 
out that TD's can support MHD surface waves (e.g. Edwin and Roberts . 1982; 
lonson , 1978; Roberts , 1981a, b; Wentzel , 1979). These waves have pro- 
perties which may make them capable of heating the corona. They resemble 
Alfven waves. Thus, to the extent that Alfven waves are capable of supply- 
ing the required energies to the corona, we surmise that surface waves 
will be capable of doing so also. But, unlike the Alfven mode, the 
surface waves are intrinsically compressive, even In a linearized small ampli- 
tude analysis. As pointed out by Kollweg (1981a), this means that the surface 
waves can be dissipated by viscosity, heat conduction, and radiation, or by 
Landau and transit-time damping in a collisionless situation. 


In this paper we evaluate the extent to irtilch surface waves can be 
collisinnally dissipated in the corona» l.e* we evaluate their dissipation 
by viscosityt heat conduction and radiation. (Collisional electrical 
resistivity can safely be neglected in the corona.) Since ve are consider- 
ing collisional dissipation, we will confine our attention to dense coronal 
regions, where collisions are frequent; coronal holes will be excluded 
from our analysis. For the simple cases considered, we conclude that the 
surface waves dissipate in a resonable distance only if their periods are 
less than a few tens of seconds, and only if the magnetic field strength 
is less than ^oout 10 Gauss. The collisional damping of surface waves may 
conceivably serve to heat quiet coronal regions, if the wave periods are 
short enough. But the large magnetic field strengths which are presumed 
to exist in coronal active region loops imply that the collisional dissi- 
pation of surface waves is too weak to heat those regions. 

It should be noted that we are not considering surface wave dissipation 
by resonance absorption (lonson, 1978; Lee , 1980; Ral and Roberts , 1982; 

Lee, Rae, and Roberts , 1982). Resonance absorption may play a role in 
heating the corona, but we will formally exclude it from our present analysis 
by considering only true tangential discontinuities, of zero thickness. If 
resonance absorption does occur, it will represent an extra dissipation 
mechanism, in addition to the collisional mechanisms being evaluated here. 
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II. Basic Egtiatlons 

The dissipation of surface waves In the corona Is analysed subject 
to a number of simplifying assumptions: 

1. Gravity Is Ignored. 

2. The magnetic field pressure In the corona Is assumed to dominate 
the thermal pressure. We will therefore Ignore the effects of thermal 
pressure on the lowest-order dynamic properties of the surface waves, i.e. 
we employ the cold plasma approximation to calculate the properties of the 
surface waves in the absence of dissipation. 

3. The surface waves are assumed to be weakly dissipative In the 

sense of k^/k^ « 1, where k^ is the Imaginary part of the wave number 
and k^ Is the real part of the wavenumber component along the background 
magnetic field, B^. This permits an approximate procedure for calculating 
k^, as follows: The wave properties are first calculated in the absence 

of dissipation. Then, given those wave properties it is possible to cal- 
culate the volumetric rate at which the waves lose energy to radiation, 
and to neat via viscosity and heat conduction. Finally, the damping length 
(or k^) is determined by equating the volumetric energy loss rate to the 
divergence of the surface wave Poynting flux. This procedure obviates the 
need to calculate a full wave dispersion relation in the presence of 
viscosity, heat conduction, and radiation, but it is limited to the weak 
damping limit. Details of the procedure are given in Section V. 

4. The surface wave properties are calculated via small-amplitude 
linearized theory. 

5. Consistent with the cold plasma approximation, the background 
magnetic field Is taken to be uniform, and In the z-direction. The magnetic 


field direction does not change across a TD (c.f. Hollweg , 1982b). 


6. The background density is assumed to vary only in the x-direction. 

Thus we will be considering planar discontinuities. Moreover, we shall 
specialize to two specific cases: i. a single TO which separates two regions 
of differing density; ii. two parallel TO's which enclose a region where 
the density differs from that external to the TD's. 

7. The background flow velocity is assumed to be zero. 

Our analysis is similar to that of Edwin and Roberts (1982). Hollweg 
(1982b). Roberts (1981a. b). and Wentzel (1979). For clarity and complete- 
ness, some of our analysis of the wave properties in the absence of dissipation 
will overlap these previous works, but we consider some new features as well. 

In the absence of dissipation, the linearized MHD induction equation 
is 

11^ - 7 X («v X B^) (1) 

where B and v denote magnetic field and velocity, the prefix *5' denotes a 
fluctuate ig wave quantity, and the subscript 'o' denotes the temporally 
steady background. The linearized momentum equation 


P 


36y 
o 3t 


- 76P + 


°oz 

4n 


3dB 

3z 


( 2 ) 


where p denotes mass density, and 6P is the magnetic pressure fluctuation. 

6P - • 6B/4n (3) 

Equation (2) has utilized the fact that is constant and in the z- direction. 
Equations (1) and (3) combine to yield 
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i »■ ' 


We will find that 6v, " 0 for Che problem at hand, and Chua (4a) becomes 

z 


3iP 

3t 



7* 6y 


(4b) 


We next differentiate (2) with respect to time. The quantities 
5B and 5P can then be eliminated with the aid of (1) and (4b). In 
virtue of assumption (6), the resulting differential equation for 6v has 
coefficients which are Independent of y, z and c. We Fourier analyse In 
those variables, with solutions of the following form: 

5Vx • (Ik^z-loJt) (5a) 

6Vy - -«Vj^yCx) sin (k^y) exp (Ik^z-lwt) (5b) 

6v^ - ®*P (Ik^z-lut) (5c) 


where k^, k^, and u are constants. If these equations are Inserted Into 
Che differential equation for 6v, we obtain 


€ «V. - k B 2 - — S2. 

kx y o 3x 


32{v 




kx 


( 6 ) 



y 3x 


(7) 


6V 


kz 


- 0 


( 8 ) 


The following quantities have been defined: 
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c(x) - 4irp^(x)u^ - 

q2 - _ e/B^2 

■ k^ - «2/v^2 

k^ - k 2 + k 2 

y * 

V.2 - B 2/Axp 

A O O 


(9) 

(10a) 

(10b) 

( 11 ) 

( 12 ) 


Equations (6) and (7) combine into the fundamental equation of this 
paper : 



dx 


)- 


c6V. - 0 
kx 


(13) 


Equation (13) is a special case of similar equations which have appeared 
previously (Edwin and Roberts , 1982; Hollweg , 7.982b; lonson, 1978; Roberts 
1981a, b; Wentzel , 1979). 

Up to now Pg(x) is an arbitrary function. We now specialize to cases 

where p (x) is constant except at the TD's, where p^(x) jumps dlscontlnuously. 
o o 

At a discontinuity, equation (13) must be solved subject to the following 
jump conditions: 


8 V, • constant 
kx 


e 36V,oj 

— *• — r “ const .mt 

q2 3x 


(14a) 

(14b) 


Equation (14a) permits both sides of the discontinuity to move together, while 
(14b) can be obtained by Integrating (13) across the discontinuity. 
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III. The Single TD 

Consider the case idiere p ■ p^, for x<0, end P^*P^. for x>0, where 

O 01 O 02 

p . end p » ere constants. We seek surface wave solutions of the form 
'^ol ^o2 

- A exp (+|qjx), x<0 (15a) 

- A exp (-Ujlx), x>0 (15b) 


The subscripts and *2* denote x<0 and x>0, respectively. It has been 
assumed that q.^>0 and (see below). The quantity A is an arbitrary 

constant. Inserting equations (IS) into (14b) yields the dispersion relation: 



(16) 


squaring both sides of (16) and using equations (10) yields 


(c-c) [k*B^(c+€)-ee1«0 

1 2 l^y o 1 2 1 2j 


(17) 


Consider first the root e » c . This requires that either p , p or 

12 01 02 


0 . 


The former case eliminates the discontinuity and thus the surface 


wave, while the latter case merely corresponds to a new steady state. This 
then leaves the following for the surface wave dispersion relation: 


k 2 b 2 

^ E +€ 

1 : 


(18) 


Defining new variables X = ^ 2 - equation (18) takes 

the form 


- . ( X-x lifl-x) 

(6+1)Y-2X 


( 19 ) 
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where 3 s p /p , in (19) we will regard 3 and Z aa given, ao that Y can be 
02 01 

determined as a function of X. Without loaa of generality wo asMiae CK0<^, 

and we confine our attention to the first quadrant of the X-Y plane. Froei (19) 

it le easily shown that Z is positive only above Che line Y " X/3, or between 

the lines Y ■ X and Y “ 2X/(l'f0); since Z is positive definite by assvsBption, 

allowed solutions are restricted to tnese two regions. However, valid surface 

wave solutions require q ^>0 and q ^>0. By combining equations (10) and (19) , 

I “ 2 — 

it can be shown that these conditions are satisfied only below the line 
Y " 2X/(l43)t and thus valid surface wave solutions arc confined to the second 
of the two regions specified above. (It is easily shown that and have 
opposite signs in this region, as required by (16). Thus restricting our 
attention to this region eliminates the extraneous root idilch was introduced 
when (16) was squared.) 

By examining the discriminant (e.g. Thomas . 1960, Section 9-‘10) of (19), 
it is easily shown that (19) generates hyperbolas. Only the branch lying 
between the lines Y - X and Y • 2X/(14B) need be considered, for Che reasons 
given in the preceding paragraph. This bra...-h has the asymptote 

Y X+Z (20) 

and it passes through the origin with a slope 2/(lHB). The phase velocity, 
w/k^, is therefore [2/(l+6)3**v^ for low-frequency waves satisfying 
u^<<Z/(l46), and asymptotically approaches v^^^ at high-frequencies. 


IV. Th« Slab of PU«m> 


Consider two TD'a located at x» >a and x* Let " p^^ for |x|>a» 
and p^ • p ^2 Uh** where p^^ and conatanta. Tha aubacripts *1* 

and *2* will refer to |x|>a and |xj<a, respectively. We consider only solu> 
tiona lAlch are localized to the vicinity of the TD'a. We therefore aamaM 
and look for solutions to (13) of the fora 


«Vj^ - Aj «cp(-jqjjx), x>a 

(21a) 

- ±Aj exp (+jqjlx), x<-a 

(21b) 


Here Aj is an arbitrary constant,* and the ± in (21b) allows for an arbitrary 
choice of parity. 

The quantity ^2^ °^w be positive or negative. If q^^>0, equation 
(13) has the folloidng two solutions in region 2 which vt distinguish with 
the labels 'Case I* and 'Case II*: 


^2 k2|x* Case I (22) 

^2 l<l2l*’ ^*** 

On Che ocher hand, if qj^<0, equation (13) has solutions (Cases 111 and IV). 

- Aj sin Ujlx, Case III (24) 

«Vj^2 • ^2 I**:!*' 


Here Aj is a constant which in ;sch case is related to A^ by applying (14a) 
at X * a or at X * -a. (Note that Cases 1 and III have odd parity, so that 


the minus sign must be chosen in (21b). Cases II and IV require the plus 
sign in (21b).) 

The dispersion relation is obtained by applying (14b) at x « a or at 
X " -a. For the four cases defined by equations (22) - (25), we have the 
following dispersion relations: 


Case I: C - -coth (|q 2 |a) 

(26) 

Case II: % ” -tanh (|q 2 |a) 

(27) 

Case III: 5 •* cot (|q 2 |n) 

(28) 

Case IV: C * -tan (jq^la) 

(29) 

c = !i iM 
^■^2 (SI 



Note that equations (26) and (27) reduce to the dispersion relation for a 
single discontinuity (equation 16) in the limit a-^ • , which is expected since 
the two id's do not interact with one another in that limit. 

Figures la - d plot vs. a^ for k^a = 2 and Pqj/Pq 2 * 

The figures are generated as follows: The quantities k^a and Pqi/Pq 2 

first chosen. The quantity [qjl* then allowed to vary, so that each value 
of Iq^ja yields one point in Figures 1 a-d. For each value of Iq^ja, the quan- 
tities e a^/B ^ and C can be calculated via equations (10a) and (26) - (29), 

2 o 

respectively. From the definition of it is then possible to determine the 
value of (e^a^/B^^) (|q^|a) Writing e^a^/B^^ in terms of Iq^j^a^ via 
equation (10a) then yields a quadratic equation for |q^Pa^. We solve this 
equation for |qj|a, and e^a^/B^^ is then determined vla(lOa). Finally, 
equation (9) is then used to determine 
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0 9 P • 1 

“!s| . (1. * 

V* "02 V 




B 


( 30 ) 


and 


2^2 




(i)^a^ 




B 


(31) 


Several features of the solutions are worth noting: 

1. For Cases I and II, we obtain (i»^>0 and k *>0 for |q 2 l<C‘)k , If P_,<(^)P 


2. For Cases I and II, and k^^-K) as Jq^Kk^. In this limit, 

it is readily shown that 


02 


(t) 


* coth(k^.)] 


,2 % 


(32) 


for Case I. (For Case II, the coth in (32) must be replaced with tanh.) 
3. For Case I with u^<" and k^^*as |q, In this 

limit 


^ :^;k2a2+kV 

* y 


(33) 


(see equations (31) and (10a). 

4. For Case II with p , < p o> ^<^<1 k^ approach maximum values as 

0 1 0 ^ z 

|q, |->^0. This maximum value is evident in Figure lb. It is readily shown 


that these maximum values are 


ci)2 « k ^v. ^ 
max y A2 




1 - <»ol''>o2> 


(34) 


2 (»o i/%2> (y)^ [n-(l-f4k^'^.~^)‘^']/2 


z , max y 


1 - (Pol/Po2> 


( 35 ) 


5. Cases III and IV exist only when < Pqj. Then 

the region between the two TD's can act as a waveguide for the fast MHD 
node, with total internal reflection occuring at the TD's. Strictly speakings 
Cases III and IV are not true surface waves > but trapped body %raves (Edwin 
and Roberts, 1982). We include Cases III and IV for coapleteness. 

6. The point indicated by the arrow in Figure Id (Case IV) is for 

Id.. 1^0. It is identical to the maximum values of u and k found for Case II 
(equations 34 and 35). 

7. For Cases III and IV, equations (28) and (29) show that 

for certain values of Iqjia* For these values of |q 2 i&> it is easily shown 


that 


(i)^a^ 


Iq^lV 


(36) 


and 


2 - 2 , 


k ‘a 
z 




k 2.2 

y 


(37) 


Some of these points are indicated by the x's in Figures Ic and Id (there 

are some additional points with k ^<0, in virtue of equation 37 , which are 

z 

not shown in the figures). 

8. For Cases III and IV, equations (28) and (29) show that for 

certain values of [qjl^. For these values of [qjla, it is easily shown that 

and k 2.*.<» with 
z ’ 


A2 



( 38 ) 


asymptotically (see equations 31 and 10a) . 
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V. Colllslonal Dlsaipation . 

As discussed In the Introduction, we wish to calculate the dissipation 
of the surface waves by viscosity, heat conduction, and radiation. We will 
do the calculation for the plasma slab (see Section IV); the case of a 
single discontinuity can then be obtained by letting a 

The volumetric heating rate due to viscosity, Is given by 

Braglnskll (1965, equation 2.28): 

Q^ls - ^ (V«6 v)2 (39) 

where It Is assumed that 6 v ~ 0 and that the Ions imdergo many cyclotron 

z 

orbits between collisions. The viscosity is due mainly to the protons. 

From Braglnskli's equations (2.U) and (2.22) we have 

n - 1.0 X 10“^® T 
o o 


in c*g*s* units; T is temperature (degrees Kelvin) and the Coulomb logarithm is 
taken to be 22 . 

The heat conduction flux is carried matnlv by the electrons. If the 
electrons are magnetlzedt the volumetric heating rate due to heat conduction 
( Braginskii ^ 1965, p.303) 


^ther 


K 



(40) 


where < is the electron heat conductivity along B . From Braginskii's 
(l,e •'O 

equations (2.12) and (2.5e) we have 


-7 5/2 

K - 8.4 X 10 T 
n ,e 0 


(41) 


in c.g.s. units; the Coulomb logaritlim is again taken to be 22. If the wave 


frequency Is very small, then the waves induce nearly adiabatic compressions 
in the plasma, and we have 


|6t| ■ (y~ 1) T^“~^|v*6y| (A2a) 

and equation (AO) becomes 

%her “ •'lue T^(V-1)2 (V-6v)2 (A2b) 

(Y Is the usual ratio of specific heats, which we shall take to be 5/3). On 
the other hand, heat conduction can smear out the temperature fluctuations of 
higher-. requency waves. In the limit of very high frequency, one has 

|6tI - n KT k |7*5vl (A3a) 

and 

Q., • n k (7»5v)^ (A3b) 

^ther o o 2 ii,e - 

where K is Boltrmann's constant. At Intermediate frequencies neither (42) nor 
(43) is applicable. We will not do the calculation at this level of detail, 
however. Instead, we shall use equations (42) if equations (43) 

if where is the angular frequency at which (42b) and (43b) are 

equal, l.e. 

( 4 ) m 

cond n K 

o 

Finally, we assume that the plasma is optically thin, and radiates 
An^ erg where n is the electron number density and A is a tempera- 
ture-dependent factor which we shall take to be lO'^^^erg cm^ at coronal 
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temperatures (e.g. Rosner, Tucker end Vslsns, 1978, Figure 10). For the 
purposes of this analysis, we consider only the radiation Induced by the 
waves. The volumetric rate at which radiation extracts energy from the 
waves Is given by 


Qrad 


But mass conservation implies 


l«n| 


no V»5v 

(I) 


(44a) 


<44b) 


and thus 

Q , ■ An ^u”^(V»6v)^ (44c) 

^rad o 

Equations (39), (42b) or (43b), and (44c) are the desired expressions 
for the volumetric rates at which the %raves lose energy to the plaana and 
radiation field. To evaluate w..ese expressions we need to know 7»4y. From 
equations (5)-(8) we obtain 

V»6v - - (e/Bjj^) ky"^ 5Vj^y cos(kyy)exp(lk^z-lwt) (45) 

The next step In the calculation requires knowledge of the Poyntlng 
flux carried by the wave. If the wave electric field Is given by 
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4E ■ - fiv X B^/c (46) 

(c Is the speed of light) end if 4v • 0, then the s-component of the Poynting 

z 

flux is 

S ■ ^oz 6v*6B (47) 

* 4tt 

But equation (1) gives 


6B ■ - k B 6v /u 

X Z OZ X 


(48a) 


and 


6B ■ - k B 6v /u) 
y z OZ y 


Inserting (48) into (47) gives 

®oz^ Lj 


(48b) 


S - ^ (5v % «v S 

z Alt (i> X y 


(49) 


Equation (49) is the desired result for the Poynting flux; it is easily shown 
that the time averages of and vanish if the solution has the form of 

equations (5). 

The damping length of the waves is calculated by equating the volumetric 
heating and radiation rates to the negative of the divergence of the Poynting 
flux. However, we first average all quantities over time and over the y- 
coordinate. And, because the waves vary in x, it is necessary to first in- 
tegrate over the x-coordinate. Thus, if we denote the e-folding distance 

for the Poynting flux by L , we have 

z 

<S_>dx 

9D * 


,%is * ’ther * <!r.d ’ 
where the ovcrbar denotes the average over time and the angle brackets denote 

the average over the y-coordinate. 
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Evaluation of the Integrals in (50) is somewhat tedious » and we merely 
quote the result. For cases I- IV we find 


L 


z 


®2.®^ N 

4n 0 ) D 


For Case I, N and D are 


(51) 


Nf ■ |qj"^(l + ky2/q^2) sinh^lq^ja 

+ (sinh 2|q |a - 2s|q Ia)/(2|q |) 

2 2 2 

+ k ^(sinh 2|q |a + 2slq la)/(2|q | 3) (52) 

y 2 2 2 

D - e 2d (sinh 2|q |a + 2slq |a)/(2|q P) 

^ 2 2 2 2 2 

+ e^2d^ Sinh2(|qja)/|qj3 


where a 


+ 1 . 


If «<ai . 
eond. 



the quantity d la 


(54) 


“^'“cond. replaced by on the right-hand side of (54). 

The damping length for Case II is given by a similar expression, which we will 

not write down; it is obtained from (52) and (53) by taking s “ — 1 and by replacing 
slnh2jq^|a with cosh2|q^|a. For Case III we find 


*^III " sln2|q^|a 

+ (-u*sin 2|q^|a + 2 |q^ |a)/ (2|q J ) 

+ k^2(„.gi„ 2jq^|a + 2lq^|a)/(2|q^P) (55) 

°III " ^*‘*2^® * 2|q^|a)/(2|q^|3) 

+ e^2j^ sln2(|qJa)/|qJ3 (56) 


where u - + 1. The damping length for Case IV is given by a similar expression; 
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it la obtained from (55) and (56) by taking u ■ >1 and by replacing aln^jq |a 

2 

with coa^la la. 

2 

The behavior of L aa a function of wave period la illuatrated In Flgurea 

m 

2a-d for Caaea I-IV. We have taken p /p ■ 0.25, k a ■ 1.0, B ■■ 10 Gauaa, 

01 02 y o 

a ■ 10®cm, n ^2 ■ 5 x lO^cra”^, and ■ 2 x 10®K. Except perhapa for 

the magnetic field atrength, theae parametera might correapond to a aolar active 
region or to a denae quiet region. The effect of radiation tuma out to be 
totally negligible for thia choice of parametera, and the damping la almoat 
entirely due to viacoalty and heat conduction. 

Figure 2 ahowa that the damping length tenda to become very large for 
long wave periods. This Is particularly evident In Figure 2a, where It Is 

seen that L « b)'** at the long periods. (This dependence on u results from 

z 

the frequency dependence of the *e 2 and c 2 factors appearing in equations (53) 

2 1 

and (56). These factors reflect the frequency dependence of the compressive- 
ness of the waves, l.e. equation (45): lo%ir-frequency surface waves are only 

weakly compressive, and thus only weakly dissipative. This fact severely 
restricts the ability of surface waves to heat the corona. A survey of para- 
meters of Interest to the corona has revealed that only waves with periods less 
that a few tens of seconds will be able to dissipate their energy before pro- 
pagating a few tenths of a solar radius. 

By the same token, equation (45) shows that the compressiveness is 
strongly reduced if B^ Is increased. We can therefore anticipate the result, 

Co be demonstrated below, that the surface wave damping becomes negllble If 

B exceeds 10 Gauss or so. 
o 

Fip^ures 2c and 2d show that the trapped fast waves do dissipate rapidly. 

If fast waves can be i^enerated in the corona, they can heat the corona. But it 
must again he recalled that the sun may not generate waves of such high frequency. 
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How do different choices of parameters affect Figure 2? A survey of 
parameters Ir the vicinity of those used In Figure 2 has revealed the following; 

1. If only Is changed: The curves In Figures 2a and 2b move up 

roughly In proportion to to the right roughly In proportion to 

ky"^. On the other hand, the curves In Figures 2c and 2d are not strongly 

modified by changes In k^, except that smaller values of k^ allow the curves 

to extend to longer periods and larger values of L (and vice-versa for larger 

2 

values of k^). Overall, we have found that values of k^a between about 0.5 

and 2 lead to the "best" surface wave dissipation In Cases I and II. For 

smaller values of k a the curves are moved up to very large values of L , 

y z 

while for larger values of k^a the curves are moved to the left to very high 
frequencies. 

2. If only a Is changed while k^ and all other quantities are held 

fixed: Figure 2a la essentially unchanged, while Figure 2b Is changed mainly 

by extending to higher (lower) frequencies If a Is Increased (decreased) . 

The curves In Figures 2c and 2d move to the right roughly In proportion to a, 
and at high frequencies they move up roughly In proportion to a^; the maximum 
levels reti^.hed by the curves in Figures 2c and 2d do not respond strongly 
to changes In a. 

3. Increasing (decreasing) B increases (decreases) L . For Cases I 

o z 

and II at low frequencies, where L « oa”**, we find that L scales as B ^ 

2 2 0 

at a given value of u). (This result Is in part dependent on our being In a 
parameter regime where equations (42) apply at low frequency, and heat con- 

duction dominates viscosity in equation (54). If T ■ 2 x 10®K, heat con- 

o 

duction will dominate viscosity if o)/k < 1500 km s"^.) At high frequencies 

in Figures 2a-d, we find that L_ « B . This means that the waves dissipate 

to 

much less efficiently in regions of strong magnetic field, and vice-versa. 


The ability of surface waves to colllslonally daop and heat the corona will 
therefore be limited to weak magnetic fields, as in quiet coronal regions. 

Our survey of parameters has revealed that the damping lengths become un- 
acceptably large If > 10 Gauss. Coronal active regions are presusMd to 
have field strengths well in excess of 10 Gauss, and this Mchanism fails 
in those regions. 

4. We have also investigated the consequences of changing n^^ 

while n Jn and all other quantities are held constant: The effects of 

changing the density are difficult to categorise in detail. But we have 

found that in general denser regions have shorter dissipation lengths, and 

are more prone to heating by colllslonal wave dissipation. 

Finally, it is useful to note that in most cases where the wave dissi- 

* 

patlon is significant, the bulk of the heating occurs in the denser region. 

To Illustrate this point, we plot in Figures 3a-d the ratio of the Integrated 
heating rate in region 2 to that in region 1, for Cases I-IV, l.e. we plot 
the quantity 

R = f o^^vis ^ther^*^* v-7) 

I ^^^is + ^ther^*^* 

vs. wave period. (For Case I, R is formally equivalent to the ratio of the 

first term on the right-hand side of (53) to the second term with A taken to 

be zero in d and d . For Case III, R is the ratio of the first to the second 
2 1 

term in (56), and similarly for Cases II and IV.) it is clear from the 
Figures that the heating occurs primarily in region 2. 
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VI Validity . 

Equation (39) for tha viacoua haatii^s rata la valid for a colliaional 
plasiM. In tha praaant contaxt» thia raquiraa 

«T <1 (58) 

P 

whara t la tha proton collision tlna. From Braginakll'a aquation (2.5i) va 
P 

find 

T -0.75 a (59) 

P P P 

for a Ck>ulomb logarithm of 22 (tha aubacrlpta 'p* danota protons). 

Similarly, equation (40) for the heating due to heat conduction la 
valid only if the electron mean-free-path, is short enough so that 

k X < 1 (60) 

z e 

If X - T (KT /m )**, where K is Boltzmann's constant, m^ the electron mass, 
e e ' e e a 

and T the electron collision time (Braginskli's equation 2.5e). we have 
e 

X - 4.9 X 103 T 2 n cm (61) 

e e e 

Figure 3 shows that the dissipation occurs mainly in region 2, and we 
therefore evaluate (58)-(61) in that region* We also note that (60- (61) will 
be more Important at lower frequencies, where the heat conduction dominates 
the viscosity, while (58)-(59) will be more Important at higher frequencies, 
where viscosity dominates heat conduction. 

Overall, (58)- (59) turns out to be the most relevant constraints for 
parameters in the vicinity of those used in Figures 2 and 3* They restrict 
the validity of our calculations to wave periods greater than 3 seconds or so* 
The left-reo.Ht portions of KlKurcs 2 and 3 siiould therefore be disregarded. 


VII. Diicus>lon 


Umvc theorlas of coronal haaclng faca tha dual ctiallanga of idantlfylng 
a wave mode %rhich la aimultaneoualy able to propagate the required energy 
flux into the corona, and to diaalpate that energy there. Hollveg (1981a) 
has suggested that MHD surface waves can In principle satisfy these require- 
nents. Their propagation is slailar to that of Alfven waves, which appear 
capable of propagating energy into the corona (Hollweg , 1981b). And they 
are compressive even in linearized theory, and therefore capable of heating 
the plasna via ordinary viscosity and heat conduction. 

In the preceding Sections we bcve discussed coronal surface waves with 
the Intent of assessing whether they can in fact dissipate rapidly enoigh 
to heat the corona. We have considered only planar geometries and two cases: 
a single discontinuity separating two regions of different density, and ttio 
parallel discontinuities which enclose a denser region. (The case of two 
parallel discontinuities admits solutions which are really trapped fast waves, 
rather than surface waves. We have Included these trapped fast wave solutions 
for completeness.) We have assumed that the discontinuities are structureless, 
so that resonance absorption need not be considered. And we have assumed that 
the magnetic pressure dominates the thermal pressure, so that the wave dynam- 
ics can be approximately calculated by assuming the plasma to be cold. (Some 
of our numerical examples have stretched this latter assumption somewhat. For 
example, the parameters in Figures 2 and 3 imply that the total (electron plus 
proton) thermal pressure is two thirds of the magnetic pressure.) 

We have carried out the analysis under the assumption that the wave 
periods are long enough for the plasma to be collision-dominated. For the 
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p«ru«t«r« of IntorMt, this roatricts tha analyala to tmva parioda lont«r 
than a faw aaconda. Wa hava alao only approxlaatad tha functional bahavior 
of tha taaparatura fluctuationa In tha praaanca of haat conduction. Our 
approxiaation ia good in tha high and low fraquancy llaiita, but wa baliava 
it to ba ganarally adaquata for tha praaant purpoaaa. 

Wa hava found that tha aurfaca wava diaaipation langtha ara aaallar 

than a faw tantha of a aolar radiua only if tha wava parioda ara laaa than 

a faw tana of aaconda. Longar pariod aurfaca wavaa hava vary long daatping 

langtha, and tha wavaa ara than not auitabla for coronal haating. Tha trappad 

o 

faat wavaa can hava abort damping langtha, of tha ordar of 10 cm or laaa, and 
thay can haat tha corona. But tha trappad faat wavaa axiat only for short 
wava parioda, 10 aaconda or laaa in our axampla, and it ia not claar if tha 
aun ganerataa wich wavaa. 

The requiremant of auch high-fraquancy wavaa ia a fundamental difficulty. 
Thera ia no evidence that tha convection zona ganerataa motiona with auch 
hl« t frequanciaa. Evan if motiona in the required frequency range ware 
generated by the convection zone, thay wouU probably diaaipata in tha photo- 
aphera and chromoaphare, and not penetrate into tha corona aa high-frequency 
waves. On tha other hand, Bruackner (1980) has reported ebsarvations of rapid 
changes in the chromosphere-corona transition region, occurring on time scales 
shorter than 20a; but it la not clear whether these motions ara waves; they 
seem rather to be Impulsive in nature. It is alao possible that high-fre- 
quency waves can ba excited on coronal loops, by the resonances pointed out 
by Hollweg (1981b), lonaon (1982 ), and Zugzda and Locana (1982); as yet there 


is no observational evidence that these resonances are excited. 


Another difficulty with this mechanisn is its extresie sensitivity to 
the magnetic field strength. Increasing slightly can drastically increase 
for the surface waves. And increasing means that the trapped fast %«ves 
can exist only at even higher frequencies, which s»kes their presence in the 
corona even less likely. On the whole, we find that the waves considered in 
this paper do not dissipate effectively if B^> 10 Gauss. Larger field 
strengths are presumed to exist in active region loops, and the mechanlaa 
investigated here falls in those regions. 

Unless the sun generates coronal surface waves with periods slx>rter than 
a few tens of seconds, we conclude that the collislonal dissipation of surface 
waves is not an effective means -of heating the corona, contrary to the suggestion 
of Hollweg (1981a). 
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Figure Captions 


Fig. 1. The wave dispersion relations, «*a*/v^ 2 * vs. k^*a^, for surface 
waves (Cases I and II) and trapped fast waves (Cases III and IV). The 
graphs are for the case of two parallel TD's, with Pqj/Pqj* 

Fig. 2. The dissipation length, L , vs. wave period for surface waves 

(Ci\ses I and II) and trapped fast wives (Cases III and IV). The graphs 

are for the case of two parallel TD's. The paraneters for all four cases are 

.4 

given In (a). For the surface waves. It Is found that L at periods 

z 

longer than a few tens of seconds. Surface waves with periods longer than 
a few tens of seconds can not h^t the corona. 

Fig. 3. The ratio of the Integrated wave energy loss In region 2 to the 
Integrated energy loss In region 1, for the case of two parallel TD's. All 
parameters are the same as In Figure 2. There Is a strong tendency for most 
of the wave energy loss to occur In the denser region. 
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